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, ^ , We show that the correlation between morphological fluctuations of two interfaces confining a 

' quantum well strongly suppresses a contribution of interface disorder to inhomogeneous line width 

of excitons. We also demonstrate that only taking into account these correlations one can explain all 
the variety of experimental data on the dependence of the line width upon thickness of the quantum 
well. 
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C3 ■ I. INTRODUCTION 

Absorption and luminescence exciton spectroscopy are one of the most important tools for studying quantum wells 
^ ' (QW) as well as other semiconductor heterostructures. Therefore, one of the most fundamental problems in physics 
of these systems is establishing connections between spectral line shapes and microscopic properties of the respective 
structures. A great deal of efforts was devoted to this problem over the last half of century, and it has been established, 
among other things, that at low temperatures the exciton line width in absorption and photoluminescence spectra in 
. quantum wells is predominantly inhomogeneousi. The shape of the spectra in this case is determined by various types 
' of disorder present in a structure, and it is currently generally accepted that the spectral width in QW is directly 
^ related to the quality of interfaces, so that the luminescence spectra provide a quick and simple quality assurance tool 
Q for QW growth^. 

O However, as it will be seen from the subsequent analysis, in spite of all the efforts, the existing theories of inho- 

mogeneous broadening of excitons in QW still cannot satisfactory explain all the diverse experimental data collected 
in this area. While it is quite plausible that the present theories suffer from several drawbacks, it is the goal of this 
^ , paper to point out at one of the important circumstances completely overlooked in the previous studies. We will 
C3 ' show here that it is not possible to explain experimental results collected by different groups for different material 
systems without taking into account so called inter-wall correlations, which describe the fact that properties of two 
\l interfaces forming a quantum well are correlated, especially for narrower wells. While direct morphological analysis 
', done with cross-sectional scanning tunnel microscopy"^, scattering ellipsometrjii and X-ray reflection measurements 

■ give convincing experimental evidences of such correlations, their role in optical spectra of QW was never studied 
' before. Another relevant example of such correlations is vertical stacking of quantum dots^, where vertical correlation 

"p^ . length is observed up to the 80 monolayer thickness. 

■ One of the main objectives of this paper is to show that the inter-wall correlations can significantly modify optical 
, spectra, and that taking these correlations into account is necessary in order to achieve even a qualitative agreement 

I ' between theory and experimental results. One of the important practical conclusions of this work is that narrow lines 
■ O [ not always mean a good quality interfaces contrary to the popular belief, but can be the result of a line narrowing 
^ ■ effect of the inter- wall correlations. 

^ ' The paper is organized in the following way. In the next section we shall provide a brief review and critical analysis 
i;^' . of existing experimental results and relevant theories. In subsequent section we will generalize the earlier theories of 
the interface disorderiSiSi^iiSiii to include effects of the vertical correlations. In the last section we will compare the 
results of our analysis with experiments. The paper is concluded by an Appendix, in which we comment on the role 
of the lateral (in-plane) correlation length in the optical spectra of QW. 
a ' 

_ _ I 

II. COMPARISON BETWEEN CURRENT THEORIES AND EXPERIMENT: CRITICAL REVIEW 



Since a QW is generally a heterostructure formed by a binary semiconductor (AB) and a ternary disordered alloy 
(ABi-xCx), there are two types of disorder responsible for the inhomogeneous broadening. One is compositional 
disorder caused by concentration fluctuations in a ternary component of the QW as well as random diffusion across the 
interface^^-^'^. The other source of inhomogeneous broadening in QWs is associated with the roughness of the interface 
caused by the formation of monolayer islands at the interfaces resulting in local changes in the well thickness iSiSiSiifliii 
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The quality of the interfaces is very sensitive to the ambient parameters of the growth process. Depending on growth 
conditions, the atoms deposited on the surface can form "islands" of various lateral sizes with different correlation 
scales. These morphological changes manifest themselves in the shape and width of optical spectra of QW excitons. 
Since both types of disorder can be ultimately traced to local changes in concentration, an accurate distinction between 
them is not a trivial task, and was first elucidated in Ref. 

Theoretical studies of effects due to compositional and interface disorder on absorption and photoluminescence 
spectra have a long history (for review articles see, for example, Refs. lll lT^ITsI and references therein.). This problem 
can be divided into two fairly independent parts. The first deals with the derivation of the random potentials 
acting on excitons in QW's from the properties of microscopic fluctuating parameters (concentrations, well thickness, 
etc.). Its main objective is to calculate correlation functions of these potentials. The second part of the problem 
consists of calculations of characteristics of excitons subjected to these potentials and in establishing relations between 
characteristics of optical spectra (line width, shape) and properties of the potentials (r.m.s. fluctuation and correlation 
length). Both these problems were carefully studied in the past, but since the focus of the present work is on the 
former we shall discuss it in more detail. 

The main object of our discussion is the r.m.s value of the random potentials, W, defined as = V (^e//(R)^), and 
its dependence on the microscopic parameters of the QW. Here Vef / is the effective potential acting on excitons due 
to both compositional and interface disorder. Current theories predict distinctly different properties for contributions 
to the potential from these two types of disorder. In particular, the dependence of W upon the width of the well, 
L, for a narrow well is predicted as being cx L for the interface disorder. For the compositional disorder, it depends 
on whether the QW is formed by a ternary alloy or a binary systemic. In the former case (as in Ina;Gai_a;As/GaAs) 
W oc L'^/^, while in the case when the QW is formed by a binary material, the dependence is cx VX- In the three- 
dimensional limit L —^ oo the contributions of the two types of disorder to W{L) are also different. The contribution 
of the interface disorder for large L decreases as while the role of the alloy disorder again depends upon the 

type of the structure. If QW is formed by a ternary alloy W only weakly, as decreases with the width before 

reaching a constant three dimensional limit. If, however, the alloy forms barriers, W decreases with L much faster 
as exp{—KoL)/L^, where kq is an inverse penetration length of ground state wave function in the barrier region. 
Although in both cases the function W{L) has a maximum at some intermediate values of L, the position of the 
maximum, and the shape of the function W{L) differ for the two types of disorder. 

While the qualitative picture of disorder induced broadening is understood rather well, attempts at a quantitative 
comparison of theoretical predictions with experiments face significant difficulties. In Fig. ^ we have collected exper- 
imental data for dependencies of low temperature photoluminescence exciton line widths on QW average thickness, 
L. All the data are for luajGai-icAs/GaAs heterostructures, and represent experimental results from several research 
gj.Q^pgi6ji7ji8ji9_ -y^^hile all the results show a non-monotonous dependence in accord with theoretical expectations, 
the maxima have different positions and, do not seem to have a regular dependence on concentration: the peak for 
X = 0.18 lies between peaks for x — 0.12 and x — 0.135. The maxima also have different heights and sharpness. For 
example, the full-width at half maximum (FWHM) for x — 0.09 and x — 0.11 have very smooth behaviors more char- 
acteristic for compositional disorder, while other data have rather sharp features more typical for interface disorder. 
Finally, the values of FWHM at large L, which are determined mainly by compositional disorder, are scattered over 
quite a broad range. 

Let us analyze Fig. ^ in light of present theories. Regarding the alloy contribution to the line width, we note 
that in quantum mechanical approaches to this problem.^'^^'^^'^*^ this contribution is completely determined by alloy 
concentration and QW width; the theories do not contain any unknown parameters that could be used to flt theoretical 
predictions to the experiments. Therefore, one can directly compare theoretical predictions with experimental results 
using the large L asymptote of the line width. Initial calculations for alloy induced disorder in the bulk^^-^^ and in 
the quantum wells^ were done in the adiabatic approximation, where Bohr's radius of an exciton, as, was assumed 
to be much larger than a correlation length of the effective potential, £c- If one applies that theory^^ to the case of an 
Ino.12Gao.88As/GaAs QW, the FWHM comes out to be = 0.36 meV for L = 150 A. This value is much smaller than the 
observed bulk valued. The same order of magnitude results are obtained in the semi-classical limit of the theory^^'« 
with the Gaussian shape of the exciton line width. One can reasonably arguai that the adiabatic approximation fails 
for a heavy hole and a light electron (m/j ^ rrie) both being subject to short-range energy fluctuations. Even for 
underlying "white-noise" disorder, the effective potential felt by an exciton has two different correlation lengths. A 
rather massive hole will average only the small volume around center-of-mass (COM) of the exciton, whereas a light 
electron is spread out over a much larger area of the order of a^. As a result, the hole will be much more sensitive to 
the compositional fluctuations, and its contribution to the effective disorder potential will be enhanced by the factor 
(M/me)^, where M = me + w/t- However, the line widths found using improved theories^'^'' turn out to be much 
larger than the experimental results (see also Fig.[71below). Thus, existing theories cannot produce an accurate result 
for the alloy disorder contribution to the exciton line width. 

It is less straightforward to compare the theory with experiments for interface roughness induced broadening because 
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FIG. 1: Experimental dependences of the low temperature exciton full width on half maximum (FWHM) on the QW average 
size, L. All the data are presented for Ina;Gai_a;As/GaAs QW's. The results were taken from the following references: Patane, 
et allii£, Kirby, et allpii, Bertolet, et allii^, Reithmaier, et allii^. 



it is difficult to separate contributions from two types of disorder in the regime of small and intermediate values of L. 
One could hope to identify the most important contribution by the slope of the W{L) dependence at small L, but, 
unfortunately, the accuracy of the existing data does not allow one to distinguish between L or L'^/^ dependencies. The 
manifestation of the interface disorder in optical spectra is clearer in the case of systems where growth interruption 
resulted in formation of sufficiently large monolayer islands. It was shown in Ref. [n] that these islands could be 
responsible for the observed splitting of the exciton spectra. Also, interface roughness has been studied directly by 
such experimental techniques as microphotolumincscence, cathodoluminescence, transition electron microscopy, or 
scanning tunnelling microscopyi^iifliiiiiSii?-^*-^^-^^-^''^-^''--^^ 

The diversity of experimental behavior for FWHM shown in Fig. presents a significant difficulty for existing 
theories of interface contribution to the line width^-^* even from the point of view of qualitative interpretation of 
the results. Indeed, the statistical properties of an interface are usually characterized by two length parameters: the 
thickness fluctuation, h, and the lateral (in- plane) correlation length, aj^. With some reservations they are often 
taken to be equivalent to the average height and the lateral size of the islands at the interface. The size of height 
fluctuations, /i, usually has a very restricted range of variations of one or two monolayers. In the semi-classical limit 
of the Gaussian shape line width, the broadening is proportional to the product ha±, which is really only one flexible 
parameter of the theory. With the help of this parameter it is possible to adjust the relative heights of the maxima 
but not their positions and particularly the sharpness of the peak. Even if one takes into account contributions from 
both types of disorder, it is still not possible to explain the variations between the optical spectra of different samples, 
a problem addressed in the last section of the paper. 

It is clear from the provided analysis that existing theories of inhomogeneous broadening of excitons are unable 
to quantitatively explain experimental data. We suggest in this paper that one of the reason for this failure is the 
neglect of inter-wall correlations mentioned in the Introduction. While this idea does not fix the problems caused by 
wrong estimates of alloy disorder contribution, we will show in the subsequent sections of the paper that it does allow 
to explain all variety of experimental results related to the properties of the curves in Fig. ^ in the vicinity of their 
maxima. 



III. STATISTICAL PROPERTIES OF THE INTERFACES AND EXCITON EFFECTIVE POTENTIAL 



A. A model of the interface disorder 



In order to make the problem tractable, we introduce standard simplifications assuming that both conduction and 
valence bands are non-degenerate, and that they both have an isotropic parabolic dispersion characterized by the 
masses me and rrih, respectively. Throughout the paper we use effective atomic units (a.u.), which means that all 
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FIG. 2: Characteristic length scales describing interface roughness of a QW. 



distances are measured in units of the effective Bohr radius = fi.^e//x*e^, energies in units of Ea.u. = n*e'^ /h^t^ = 
2 Ry, and masses in units of reduced electron-hole mass ^l* , where 1/^* = + l/rnf^. In this notation me,h = 

m*f^/lj,*, where m* ^ are effective masses of an electron and a heavy hole. We will choose z-axes in the direction 
of growth of the structure (vertical direction). The plane perpendicular to this direction is the lateral plane (see 
Fig. 121 . We measure electron and hole energies in the QW from the conduction and valence band edges of the barrier 
respectively. Then the potential of a QW with interface roughness is given by 

C/e,h(r) = -Ve,h [e{z + L/2 - (:e, y)) - 6{z - L/2 - 7n{x, y))] « U^^l{z) + 5Ul^lf{v), (1) 

where 9{z) is a step-function, Ve^h are differences in off-set band energies, and 

Ufliz) = -V,.u[(>{^ + Ll2)-e{z-Ll2)], (2) 
5Ulf{v) = V,Mx,y)5{z + L/2)-i^2{x,y)5{z~Ll2)]. (3) 

Random functions 771^2(2^,2/), with zero mean, characterize the deviation of the ith interface from its average position. 
The perturbation expansion of the 0-function is justified because an interface roughness is almost always small for 
typical parameters in semiconductor heterostructures. The presence of two functions, 771.2(2;, y), distinguishes Eq.Q 
and Eq.lO from the respective equations of Reft^ where the roughness of only one interface was taken into account. 

The statistical properties of the interfacial roughness in multilayered systems can be characterized by the height- 
height correlation functions 

{rl^{{PlH{{p2))^h'f.,C{\Pl-p2\l (4) 

where h is an average height of interface inhomogeneity, and (...) denotes an ensemble average. We assume here that 
the dependence of both diagonal and non-diagonal correlations on the lateral coordinates p is described by the same 
function C,{p). For diagonal elements fa = 1, and the respective functions describe lateral correlation properties of a 
given interface {self- correlation functions). Non-diagonal elements with i ^ j introduce correlations between different 
interfaces; the respective quantity /i2(i/(T||), which can be called a cross- or vertical- correlation function, is a function 
of the average width of the well and is characterized by the vertical correlation length (T|| (a subscript "||" denotes 
that the direction of the vertical correlation is parallel to the direction of growth). 

The effect of the inter-wall vertical correlations has been previously considered in studies of conductivity of thin 
metallic film322i2&. To the best of our knowledge, in all previous microscopic studies of the exciton line shape in 
optical spectra of QWs, these correlations were omitted. Such an approximation is valid for wide QWs, but in the 
case of narrow QWs the vertical correlations are experimentally confirmed^i^i^iSii and should be taken into account. 
In the limit L/a\\ <C 1 it is reasonable to assume that the inter- wall correlation function /12 tends to unity, which 
means that for the very small separation between the interfaces one random surface spatially repeats the pattern of 
the other. (See Fig. 12). As we shall see below, the effect of interface disorder in this case tends to cancel out at least 
in the first order of the perturbation theory. For the sake of concreteness we will assume below that 

/12 = exp(-LVafl). (5) 

The value of the vertical correlation length a± depends on the growth process. In the following analysis we will also 
assume the Gaussian form for the lateral correlation function: 



C(R) = exp (-i?V2^i) 



(6) 
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The limit (t± ^ corresponds to the "white-noise" correlator 

C(R) = 27rcri(5(R). 



(7) 



Following the standard procedurei^ described in numerous papers we derive the Schrodinger equation for the center 
of mass (COM) exciton motion subjected to an effective random potential, 



+ [/ ,.(R) 



with Ug^f f (R) given by 



f/e//(R)= / iSU, + dUh)(b\p)xi{ze)xi{zh)d'pdz,dz,, = U,{R) + Uh{R). 



(8) 



(9) 



where 6Ue,h are defined in Eq. Q, p = pe — Ph, and R — [mePe + rnhPh) / M. For symmetric QW, x{~L/2) = x(L/2), 
and from Eq. ||2J) we obtain 



C/e,^(R) = l^e,/.xL(V2) j [r/i(R±/3,,eP)-%(R±/3h.eP)]0'(p) rfV 



(10) 



Here [3h.e — 'm-h,e/M . Hereafter we will omit an explicit dependence on L for the electron and hole wave functions 
x1 /i(i/2) always assuming that their values are taken at the position of the interface. The correlation function for 
the effective potential Uef / (R) can then be expressed as 



where 



T ■ 



(C/(Ri)C/(R2)) = + Thh + 2T,h, 



2h'v^xt{i-fi2{L)) / dWy(p)'/''(p')C(|R-/?,(p-p')l), 



2h'VeVhxlxl{l ~ fi2{L)) / dVrfW(p)r(p')C(|R-/3/.P + /3eP'|), 



(11) 

(12) 
(13) 



where i, j — e or h and R = Ri — R2. 

In all these expressions, the terms in front of the integrals determine the dependence of the correlation function on 
QW width L, while the integrals themselves determine the spacial correlations in lateral dimensions. Each of these 
terms can be presented in a form 



T,j{L,R) = 2h'V,VjF,,{L)G,,{R), 
P.jiL) = X?X?(l-/i2(i)), 



(14) 
(15) 



G.,(i?) = I d'pd^p'<p'{p)q^^p')Cm-pJP'+^3^p\), (16) 
where indexes i,j take values e and h (note that the order of these indices in the integrand is important). 



B. Dependence on well thickness 

Let us first focus on the functions Fij (L) which determine the dependence of the total correlation function on QW 
width L. For the sake of concreteness we consider in detail function Fhh- The analysis of other functions is similar. 
The dependence of Fhh on L comes from two factors. The first factor is the fourth degree of the electron QW wave 
function, xt^ calculated at the interfaces, z = ±L/2. This dependence in Ino.12Gao.88As/GaAs QW is shown in Fig.|2| 
by a dashed line. It is easy to understand the behavior of xt for the cases of large and small width. The characteristic 
scale here is given by l/non = 1/V2m/J4^, since this scale determines the number of energy levels in a finite QW 
(iV = 1 + [KhL/n]). 

For finite QW the ground state wave function has a piece-wise form 

/ ^cos(fcz), Z<\L/2\ 
> ~ \ Bexp{-Kz), z>\L/2\, ^^^^ 
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FIG. 3: Dependence of the correlation function of the hole-hole interface roughness potential on the well width in an 
Ino.12Gao.88As/GaAs QW. The solid line is the function Fhh{L). The dashed line is the fourth degree of the hole wave function 
at the interface. The dotted-dashed line shows approximations of this function for small and large L, given by Eqs. (1191 and 
I2UI I. The maximum of Xh{L)'^ is approximately located at I/koh = l/\/2mfiVh. 



where k = ^/ Kq — k'^ and k is the ground state wave vector. From the normahzation and matching conditions, one 
can readily obtain the value of square of the wave function at the interface 



i + Kyk'^ + KLKi/{2k^) 



(18) 



In the case of a wide QW, there are many discrete levels in the well, and the well is "almost infinite": k w tt/L — 
27r/i^Ko, K « kq, and 

^2 ^ ^^'^0 ^ ± no) 

Since for large well widths the vertical inter-wall correlations are negligible, fu 0, Eq. (|19|l determines the total 
decrease of the interface correlator with increasing width. 

In the opposite case of very narrow QWs, there is only one shallow level, which can be determined from the 
(5-functional potential approximation: k ^ kq — L'^Kq/8, k KqL/2, and 

- (20) 

If one neglects vertical correlations, this expression describes the suppression of the interface disorder for narrow wells 
because of a decreased portion of the hole (electron) wave function inside the well. This result was obtained earlier 
in, for instance, Ref 22. Inter-wall correlations, however, significantly modify this dependence. For lengths smaller 
than the vertical correlation length a\\ we have 

(21) 

where 7 is determined by the form of the inter- wall correlation factor /12. For example, for Gaussian or Lorentzian 
dependences of /i2(£), the parameter 7 = 2, while for the exponential form of this function, 7 = 1. Thus we obtain 
that, in narrow QWs, interface correlations are strongly suppressed by the factor 

F,j{L)^L^+\ L<an,l/Ko. (22) 



While the transition between two asymptotic behaviors of x is determined by the parameter kq, the behavior of /12 
depends on the correlation length cr|| , which is a completely independent parameter. Experimental data suggest that 
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it is quite possible for a\\ to be much larger then Kq. In this case, inter- wall correlations can affect not only the L ^ 
asymptotic of W , but also its behavior at L ^ kq. Instead of behavior one would have a much slower decrease 
of W with L-.Woi 

From the behavior of the wave function (|18|l we can see that the maximum of Fhh{L) is reached in the vicinity of 
L ~ l/noh- However, as the previous analysis demonstrates, the vertical correlation function /12 can significantly shift 
this position; it can also change the height and shape of the peak. Thus, the presence of the inter- well correlation term 
1 ~ /12 in function F^h can naturally explain all variety of experimental results related to the properties of the curves 
in Fig. ^ in the vicinity of their maxima. The graph of function F^h with inter-wall correlations taken into account 
is shown in Fig. |2| One can see that these correlations indeed significantly affect the shape of this function. While 
we only discussed properties of hole-hole correlator, it is clear that behavior of the electron-electron and hole-electron 
terms is similar. 



C. Dependence on lateral correlation length 



In order to investigate the results of the interplay between lateral and inter-wall correlations, let us now analyze 
the lateral correlation functions G^- (i?) . Their behavior is determined by the ratio of the potential correlation length, 
cr_L to the average size of the exciton in a plane, as well as by dimensionless parameters /3e and Ph- In order to 
evaluate the respective integrals, we chose the normalized ground state function of exciton relative motion in a 
quasi-two-dimensional forr^ 



m = ^^oxpi-p/X), (23) 

where A is a variational parameter that indicates the average exciton size. 

For the ground state function, Eq. H23() , and the height-height correlation function, Eq. © , the lateral dependence 
of the correlator Gy (i?) can be presented as a function of two parameters Gy (i?) = Gy (i?; yj, a), where 

and a — min (/3i, /3j)/ max (/3i, /?j). Parameter a is equal to unity for the electron-electron and the hole- hole corre- 
lator. For the cross term Geh this parameter is equal to nie/mh, which is much less than unity for the majority of 
semiconductor materials. This fact allows for additional simplifications when evaluating the integrals. The parameter 
Ueih) defines the ratio of the renormalized lateral correlation length Pe{h)^ of the effective hole (electron) potential to 
the original correlation length of the interface fluctuations. For holes with rrih 3> rrie this renormalized correlation 
length is much smaller than the corresponding length for the electrons. The latter implies that the interface disorder 
has a bigger impact on the holes compared to the electrons. This is also true for the contribution of alloy disorder to 
the broadening^^. Thus, the lateral correlation function can be rewritten in the following form: 

G,,{R;yj,a)^^ / dVV exp [-2y,(p' -f a^V)] exp (-|R - p' + pp) . (25) 

We shall focus on its value at the origin, Gij{0;yj,a), which determines the variance of the potential, W. The 
expression for W in this case can be presented in the following form: 

= 2h^ [1 - h2{L/ai_)] [V^xiGhhiO; ye, 1) + 2V,VhxlxlGehiO; yn, a) + K'XeG'ee(0; yh, 1)] ■ (26) 
The calculation are easier for the cross-term Geh, because we can take advantage of smallness of a <C 1: 

GehiO;yh,a)^Ayl f dttexp{^t^ ~ 2yht) + 0{a^) = 2yl/{l + af - 2yleMyl)V^i^ - ^^^iiVh)) + Oia^), (27) 



where erf(y) is the error function. The hmction H27|) is shown in Fig. ^ It has the following behavior for small and 
large y: 

. ^ - 2^yl + ■■■ , y,, < 1, 
Geh(0;y^,0)« <( ^ 3 (28) 

1-^, 
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The calculations are more cumbersome for the electron-electron and the hole-hole contributions Gii(0, y^, 1). It is 
convenient perform a transformation to the new set of coordinates which reflect the symmetry of the integral^^: 



S = pi+P2, t = pi~p2, 



Pi2 = ypi+ Pi~ 2pip2 cos ( 



(29) 



After this transformation the calculations for Gu are reduced to the one-dimensional integral 



ds e-^y^" 



2 s' 



g2 erf (is) 
2i 



(30) 



The term in square brackets in the integrand of Eq. (|3U|I is a smooth function that behaves as ^ (8/3)s'^ for small s, 



and changes its behavior to 
the correlator Gm(0; yj, 1): 



2s for s > 1. These dependencies allow one to estimate the asymptotic behavior for 



GuiO;yjA) 




y, « 1 
Vj » 1- 



(31) 



The first terms in the series expansions of Eqs. (|28|1 and H31|) for small y corres pon d to the white-noise limit of the 
height-height correlation function, Eq. Q , and were obtained previously in Ref. [SJ. They can be readily derived by 
substitution of a 5-function instead of the last exponent in the integrand of Eq. H25|l . The consecutive terms in these 
formulas describe deviations from the white-noise model in the case of the short-range correlations. The dependence 
of correlators Gij{0;y, a) on the respective parameters yi is shown in Fig. 01 from which one can see that corrections 
to the white-noise approximation become significant even at relatively small values of the correlation length, a± . 

For the "white- noise" interface roughness, Eq. Q, the analytical results can be obtained for the more elaborate 
exciton ground state trial function^ 



Hp) 



2 exp(7) 
v/27rA2(l+7) 



exp 



(32) 



which more accurately takes into consideration the three-dimensional character of the exciton. The parameter 7 = 
2d/X determines the ratio of the finite average distance d between the electron and the hole in the QW to the two- 
dimensional Bohr radius A. In this case, one can obtain for the most important hole-hole correlator (a = 1) the 
following expression 



G(0,a||,{A,7},/3,l) 



o-jl 1 -h 27 
/32A2 (l-f7)2 



(33) 
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This result formally coincides with the short-range limit of Eq. (|31|l after introduction of the renormalized 

effective Bohr radius A — A(l + 7)/ -^/l -I- 27. Since parameter 7 in this expression is usually less than or of the order of 
unity, this renormalization is not significant, and we can conclude that an approximation of the exciton wave function 
by Eq. H23|l gives reasonable results at least for the short-range interface disorder. Collecting together all results for 
the variance we obtain in the limit of short-range interface disorder 



W^^2h^ [l-/i2(L/(T[|)] 



A2 



v,txt 



PI 



(34) 



Apart from the vertical-correlation factor 1— /12 this expression coincides with results obtained in Ref.jSJ- As expected, 
in the case when holes have a significantly larger mass than electrons (the typical situation for In2:Gai_a;As/GaAs 
or A^Gai-ajAs/GaAs quantum wells), the hole-hole term in square brackets of Eq. (|34() dominates. In the opposite 
limit of the long-range interface correlations the result is 



= 2h^ [I - /i2(i/a||)] {Vhxl + Vexl 



(35) 



which agrees with the conclusion of Ref. obtained for a different model of the interface disorder: in the regime 
of the long-range correlations the distribution of the effective potential reproduces the distribution of the interface 
roughness. 



IV. COMPARISON WITH EXPERIMENTAL RESULTS 



In order to compare calculations of W with experimental absorption spectra one needs to evaluate dynamics 
of excitons in a random potential with given correlation properties. This problem was intensively discussed in 
literaturei*i^ii^*2i and we are going to use the results of the cited papers in conjunction with our analysis of the 
effective potential. There are two main models of exciton dynamics in a random model. One of them treats excitons 
quantum-mechanically in the limit of negative and large energies, while describing a most important intermediate 
region using an interpolation proccdur o^'^i^^ . In this approach the absorption line has an asymmetric shape with the 
linewidth proportional to . In the second approach excitons are treated semiclassicallj*i; if the underlying compo- 
sitional or interface roughness disorder is described by the Gaussian random process, then the shape of the exciton 
line is also approximately Gaussian with FWHM equal to A = 2^/2 ln(2)VF. A transition between quasi-classical and 
quantum regimes of exciton dynamics is determined by a parameter v = W/K^- Here Kc = f? I2MI\ is a kinetic 
energy of an exciton confined in a spatial region of size ic-, where Ic is a suitably defined correlation length of the 
random potentiate (also see Appendix). The quantum limit corresponds to the case ^ 1, while the semiclassical 
approximation is valid when 3> 1. 

In order to compare results obtained here with optical experimental data, we will make use of the semiclassical 
theory of exciton absorption, which according to Ref. can be reliably applied to the situation under consideration 
(typical material parameters for an In2.Gai_a;As/GaAs QW yield v ^ h for a ^ — 2aiat, and v increases with increase 
of the correlation length). Using the semiclassical relation for A we plot the dependence of the interface roughness 
contribution to the exciton line width as a function of the well width in Figs. |S1 and Figure [S] represents curves for 
different lateral correlation lengths ("island sizes"), while Fig.^lshows how these results are modified by the inter-wall 
correlations. The first thing to notice is that the changes in the lateral correlation length affect the height but not 
the position of the maximum and the shape of the curve. Comparison with the experimental data reproduced in this 
figure shows that an increase in cr_L drives the curves away from the experimental results. This has to be compared 
with the results of incorporating the inter-wall correlations. Fig. El An increase in CT|| not only significantly reduces 
the height of the curve maximum, but also shifts its position toward larger values of L, and widens it. 

It would be interesting to try to fit the experimental data presented in Fig. ^ with the results of our calculations. 
To this purpose one also needs to know the contribution of alloy disorder to the total line width. Making use of 
the model of short-range compositional disorder in a QWi^ without the adiabatic approximation, one can obtain the 
following formula for the alloy disorder induced variance^: 



87rA2 



-55- / Xh{z)'^ dz + Sanae I Xh{zYxe{zf dz + I Xeizfdz 

Pe J-L/2 J~L/2 Ph J-L/2 



(36) 



where aiat is the lattice constant, and ae^h = dVe.h/dx characterize the rate of the shift of conduction and valence 
bands with composition, x. The formula is given for the case of a QW made of a ternary alloy. The semiclassical theory 
of the exciton line width again yields : A = 2^/2 ln(2)M^, while the interpolation procedure for the quantum limit 
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gives A « 0.59MVr^. The Ino.12Gao.88As/GaAs QW is intermediate between these two limits, since v ^ 1. In Fig.[7| 
we present both A dependencies on well thickness. Unfortunately, as we can see, there exists a strong discrepancy 
between theoretical estimates of the contribution from alloy disorder and experimental results. It is not the goal of 
this paper to uncover the causes of this discrepancy. However, common wisdom tells us that in the bulk limit of 
very wide QW (L > as) only alloy disorder contribution should survive. The simplest way to adjust the theory is to 
introduce a phenomenological scaling down of l^comp to the value that should coincide with experimental results in 
the limit of large L asymptote. Although at present the reason for such re-scaling is unknown, it is hard to imagine 
that the proper theory of alloy disorder contribution will change the dependence of the variance on well thickness, L, 
determined by the integrals in Eq. H36() . Since our main purpose is to elucidate the role of the inter- wall correlations 
rather than to revise existing theories of the alloy disorder, we carry out this operation keeping in mind its purely 
technical nature. The results of the best fit performed in this way are shown in Fig. [71 along with the best fit values 
of lateral and vertical correlation lengths. Performing a similar fitting procedure for other experimental dependencies 
of FWHM on well thickness shown in Fig. ^ we obtained the following values (normalized on lattice constant aiat) 
for lateral and inter-wall correlation lengths in In^^Gai-a^As/GaAs QWs: for x = 0.09 (T_l = 1, cr\\ = 20, for x — 0.18 
CT_L = 1, (T|| = 10, for X = 0.11 (T_L = 1, CT|| — 5, for x = 0.12 tTj_ = 1, = 3, for x — 0.135 cr_L = 4, u\\ — 13. 
The most important result of this exercise is the demonstration that the fit would not be possible at all without 
taking into account the inter-wall correlations. We would also like to stress that it is not possible to achieve a good 
agreement with experimental results by omitting the inter-wall correlations, and using only alloy disorder scaling as 
an additional fitting parameter. We conclude, therefore, that the inter-wall correlations play an important role and 
must be taken into account when interpreting experimental results. 

V. CONCLUSION 

In conclusion, in this work we address the influence of vertical inter-wall correlations between rough interfaces on 
the exciton line shape. We show that the presence of these correlations strongly suppresses the interface disorder 
contribution into inhomogeneous broadening. The latter means that for narrow quantum wells it might happen that 
the exciton line width tells more about the quality of barrier material than about the quality of interface contrary 
to what is often claimed in the experimental literature. On the other hand, the differences in inter-wall correlation 
lengths can account for the variety of positions, strengths and sharpness of FWHM dependence on the well width for 
experimental data, obtained by different research groups. 




30 60 90 120 150 

L (Angstrom) 

FIG. 5: Dependence of the interface roughness induced broadening on the perpendicular correlation length ("island size") (t_l. 
The correlation length is given in terms of a number of lattice constants (a;at = 5.869A). For all curves the composition 
concentration is a; = 0.12, and the vertical correlation length parameter is fixed by a\\ = aiat- The dotted line is the rescaled 
alloy disorder contribution. The experimental data, shown by triangles, are taken from Ref. and are presented here for 
comparison only. 
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FIG. 7: The thick sohd hne is the best fit for the total line- width displayed together with experimental results (triangles) from 
Ref. llSL The separate contributions of the interface disorder (Ai„i) and of the rescaled alloy disorder (Aanoy) are also shown. 
The parameter-free theory for alloy disorder contribution discussed in Ref. ^3 is shown by the dashed-dotted line (without 
adiabatic approximation). The semiclassical limit of the alloy disorder contribution^ is shown by the dotted line. The fitted 
parameters are a\\ = Sajat, cr± = luiat- 
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APPENDIX A: DEPENDENCE OF THE LATERAL CORRELATION FUNCTION ON R 



Calculation of the cxciton absorption line shape is equivalent in the dipole approximation to the estimation of the 
optical density (OD) function^: 



A{e) 



E 



5{£- £i) 



(Al) 



where and '>pi{R) are corresponding energy and wave function of the Schrodinger equation ||SJ| for exciton COM. 
The shape of A{e) depends on the strength of disorder. The latter can be roughly measured by parameter = W/Kc, 
where W = y/ {Uef f{RY) is a variance of the potential energy induced by fluctuations and Kc = h'^/2M£^ is the 
kinetic ("correlational") energy of the exciton. The parameter £c determines the confinement of the exciton COM 
wave function. It can be extracted from knowledgc"'^^ of the lateral dependence on R for the correlation functions 
Gij{R, Uj, a) since 



d^R'{U{Il)U{Il~R'))/W\ 



(A2) 



For the "white-noise" height-height correlator one can show that all distances are scaled by factors (A/3i)/2, i.e. 
R — 2R/{XPi). For the cross-term correlator in the limit of heavy holes and light electrons the result is simple again 



Geh{R;yh,a) w GehiR;yh,0) = 2ylexp{-R), 
while for the diagonal terms we have Gii{R;yj,0) = ?/|/2/(i?i), where 



4 

f(R) = - 



dd 



dppexp — p — y R"^ — 2pRcos 



Function f{R) has the following limits for small and large distances: 

■exp(-^), ^>1. 



fiR) 




(A3) 



(A4) 



(A5) 



Thus, even though the initial correlator of interface fluctuations was of the white noise type, the effective noise for 
the exciton potential is colored with exponential tails and with a correlation length of the order of the exciton radius. 
Similar results were obtained earlier for the bulk compositional fluctuations'^^ and for the island model of interfacial 
roughnessii. One can show that such exponential asymptotes also persist for the long-range Gaussian correlator, 
Eq. jni, when R/y » 1. Since 



Gee{R) 



Ghh I R 



(A6) 



we can see that the lateral part of the electron-electron correlation function is suppressed by the factor {me/mhY 
but has larger correlation length. Overall, the localization length, is determined however by the hole-hole term: 
ic ~ A/3e. The normalized lateral correlation functions Gij{R) are shown in Fig. |S1 
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FIG. 8: Semi-log plot of the normalized lateral correlation functions Gij{R;yj,a). The normalization is chosen in a way to 
have the hole-hole correlator equal to unity at i? = 0. Data are given for yh = 0.1. The ratio nie/mh = 0.178 was chosen to 
depict the realistic case of InGaAs/GaAs QW. The dotted-dashed lines for hole-hole correlator are limits of small and large R 
[See Eq. ^] 
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